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Comment on "Ising Spin Glasses in a Magnetic 
Field" 

In a very interesting paper Houdayer and Martin 
analyze the T = 3d EA spin glass with a magnetic field 
-B. By using a new, powerful method, they determine an 
effective critical field B c (given by the crossings of their 
function r, the rate at which ground states change) as 
a function of the lattice size L. They use their results 
to deduce that the model is behaving like in the droplet 
approach and not like the mean-field theory, that implies 
the survival of a phase transition for finite field B. We 
show here, by using some unpublished data that the 
very interesting method and numerical results of Q] are 
completely compatible with the behavior implied by the 
Replica Symmetry Breaking (RSB) theory. 

In ||] we have introduced a dynamical approach that 
has allowed us to determine the value of the minimal al- 
lowed overlap, q m in, for the 4cZ EA model in a non-zero 
magnetic field. By establishing that this value is differ- 
ent from q ave , the average overlap measured on medium 
size samples, we have shown that also in Ad the phase 
transition to a RSB phase survives after applying a mag- 
netic field. We use here analogous data for the 3d case 
P], that also clearly show the survival of the Almeida- 
Thouless line, to interpret the result of M. 
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FIG. 1. q(t A ) versus t A at T = 0.5, B = 0.4. 

In our approach we use a large lattice (for example 
L = 32), and we select an annealing schedule to run a 
Metropolis dynamics. We select a set of decreasing val- 
ues of the temperature {T„}, starting from high T 3> T c , 
and we spend n ■ 2 k updating steps at each T n value. We 
compute q(tA), where = 2 fc , for increasing k values. 
The limit for tA — * oo of q{fA) is just q m in- The value of 
q ave is computed with usual equilibrium tempering sim- 
ulations on smaller lattices. We show in figure (|l|) what 
happens in 3d at B = 0.4 and T = 0.5, where RSB is 
clear: the fit of the data is very precise, and gives a value 
of q m in that clearly differs from q ave . Ref . [jjj deals with 
T = 0, so in figure ^ we extrapolate q m in(T) at B = 0.4 
down to T = 0. A simple power fit is very good, giving 


q m in{T = 0) = 0.74. In terms of the percentage of dif- 
ferent spins x of § this means that at B = 0.4 we find 
x ~ 0.13, corresponding to a jump with overlap 0.74: 
in the conditions of Q (where x — 0.15) one expects as 
L — > oo a crossing of the r curves close to B c = 0.4. 



FIG. 2. q m in(T) and best power fit. 

This is well compatible with the data of Q : their cross- 
ing point moves from B ~ 2.1 to B ~ 1.2 when L goes 
from 3 to 8: in these conditions it is difficult to discrim- 
inate between the droplet B c = and the B c ~ 0.4 that 
follows from RSB. RSB predicts that r ~ at small 
B, which is not far from the observed behavior. A main 
feature of short range models is that one needs very large 
L values to determine the asymptotic value of q m in- the 
value of q m in becomes sizably different from zero already 
for a small B (q m i n ~ B 3 in 4d against q m i n « Bs in 
d = oo). This implies that the width at L — oo of P(q) 
becomes much smaller as soon as B is switched on, and 
finite volume effects become much more important. This 
effect is exacerbated by the presence of the very large, 
low q tail of P(q) at L not too large (see j| for the anal- 
ogous Ad case) , that disappears very slowly for increasing 
L. The new method of § is very promising, and going to 
larger lattices will probably allow one further verification 
of RSB for finite d spin glasses. 

These simulations have been run at ZIB-Berlin. We 
thank H. Stueben for his crucial and extensive help: with- 
out him this work would have been impossible. We ac- 
knowledge W. Janke and D. Johnston contribution in 
setting up this project. 
E. Marinari, G. Parisi and F. Zuliani 
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